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Abstract 

D- dimensional cosmological model describing the evolution of a multicomponent perfect 
fluid with variable barotropic parameters in n Ricci-flat spaces is investigated. The equations 
■ of motion are integrated for the case, when each component possesses an isotropic pressure 
with respect to all spaces. Exact solutions are presented in the Kasner-like form. Some 
explicit examples are given: 4-dimensional model with an isotropic accelerated expansion at 
late times and (4 + (i)-dimensional model describing a compactification of extra dimensions. 
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^ ■ 1 Introduction 

The necessity of studying multidimensional models of gravitation and cosmology [1, 2, 3] is 
motivated by several reasons. First, the main trend of modern physics is the unification of all 
known fundamental physical interactions: electromagnetic, weak, strong and gravitational 
ones. During the recent decades there has been a significant progress in unifying weak and 
electromagnetic interactions, some more modest achievements in GUT, supersymmetric, 
^ \ string and superstring theories. 

Now, theories with membranes, p-branes and more vague M- and F-theories are being 
created and studied. Having no definite successful theory of unification now, it is desirable to 
study the common features of these theories and their applications to solving basic problems 
of modern gravity and cosmology. Moreover, if we really believe in unified theories, the early 
stages of the Universe evolution and black hole physics, as unique superhigh energy regions, 
are the most proper and natural arena for them. 

Second, multidimensional gravitational models, as well as scalar-tensor theories of gravity, 
are theoretical frameworks for describing possible temporal and range variations of funda- 
mental physical constants [4, 5, 6, 7]. 

Lastly, applying multidimensional gravitational models to basic problems of modern cos- 
mology and black hole physics, we hope to find answers to such long-standing problems as : 
singular or nonsingular initial states, creation of the Universe, its flatness, creation of matter 
and its entropy, acceleration, coincidence and cosmological constant, origin of inflation and 
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specific scalar fields which may be necessary for its realization, isotropization and graceful 
exit problems, stability and nature of fundamental constants [5], possible number of extra 
dimensions, their stable compactification etc. 

The discovery of the accelerated expansion of the Universe and the fact that the flat 
Friedmann model with the cosmological constant or quintessence now fits best the set of 
different observational data, created the problems of dark matter and dark energy. This is 
a real revolution in modern physics as we do not know now what really the dark matter 
(0.20 of 0.30) and what is the dark energy (0.70) of the total energy are. Even attempts 
to explain it via the cosmological constant or quintessence seem to change one puzzle with 
another one as necessary vacuum properties or exotic scalar fields with or without strange 
potentials are still waiting to find their place in rigorous theories, not speaking about their 
real experimental confirmation. 

In [9] we showed that the cosmic acceleration and coincidence problems may be solved 
by using an x-fluid as a quintessence and a viscous fluid as a normal matter. Viscosity of 
the normal matter can be explained by its own multicomponent structure. As is known, 
a mixture of different fluids admits a description as a single viscous fluid. We adopted 
the "second equations of state" in the form of some special metric dependence of the bulk 
and shear viscosity coefficients. These "second equations of state" generalize the so called 
"linear dissipative regime" in FRW world model, when the bulk viscosity coefficient is linearly 
proportional to the Hubble parameter. We studied D-dimensional homogeneous anisotropic 
cosmology, which allows to describe the dynamical compactification of the extra dimensions 
(see, for instance, our paper [10] on viscous cosmology). Other 2-component models in many 
dimensions also having the acceleration were found: with the cosmological constant in [11], 
with a perfect fluid in [8], with 2 non-Ricci-flat spaces [12], with p-branes and static internal 
spaces in [13], with scalar fields having exponential potentials in [17] and in four dimensions 
with a perfect fluid and a scalar field with the exponential potential in [15, 16] using methods 
developed in our multidimensional approach. 

Here with the same aim the D-dimensional cosmological model describing the evolution 
of a multicomponent perfect fluid with variable barotropic parameters in n Ricci-flat spaces 
is investigated. The equations of motion are integrated for the case, when each component 
possesses an isotropic pressure with respect to all spaces. The exact solutions are presented 
in the Kasner-like form. Some examples are given. The first example is 4-dimensional model 
with the Kasner-like behavior near the initial singularity and an isotropic accelerated expan- 
sion at late times. The other example is (4 + <i)-dimensional model describing a contraction 
of the internal space accompanied by the expansion of the external space at early times. 



2 The model 



Following papers [18]- [21] we consider the metric 



n 



ds 2 = -e 2 ^dt 2 + ^exppx^ds 2 , 



(2.1) 



i=i 



on D-dimensional space-time manifold 



M = RxM 1 x...xM, 



(2.2) 
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where dsf is a metric of the Ricci-flat factor space Mj of dimension d i: 7(t) and x l (t) are 
scalar functions of the cosmic time t, = expfa;*] is the scale factor of the space Mj and the 
function 7(2) determines a time gauge. The synchronous time t s is defined by the equation 
dt s = exp[y(t)}dt. 

Under this assumption the Ricci tensor for the metric (2.1) has the following non-zero 
components 



= e-ME^) 2 + 7o-77o 

\i=l J 

K: = e-^[x l + ^( 7 o-7)]C 



with the definition 



7o = Y.d iX \ (2.3) 



i=i 



Here indices m 8 and rii for i = l,...,n run from (D — J2]=idj) to (Z)]- =i (£,-), where 
D = 1 + £? = i di = dimM. 

We consider a source of gravitational field in the form of multicomponent perfect fluid. 
In comoving coordinates the energy-momentum tensor of such a source reads 



m 

iM(s) 



t n = H T JV ' ( 2 - 4 ) 

s=l 

(< (s) ) = diag \-p^(t),p[ s \t) : . . . ,ri s) (t), . . . ,p«(f), . . . ,pW( t ) , (2.5) 



ditimes dntimes 

Furthermore, we suppose that the barotropic equation of state for the perfect fluid compo- 
nents is given by 

P S' ) (t) = (l-/iS' ) (x))pW(t), s = l,...,m, (2.6) 
where variable barotropic parameters are given by 

h^\x) = j^'Hx), i = l,...,n (2.7) 



with an arbitrary smooth function $( s )(a;) on IR n . 

The equation 
tensor (2.5) reads 



The equation of motion XJ mTq 1 ^ = for the perfect fluid component described by the 



pV + Y,d i #(p w +pi' ) ) = 0. (2.8) 
i=i 

After using equations of state (2.6) via (2.8), integrals of motion may be obtained in the 
form 

A (s) = p {s) exp J 27o _ $ ( S ) ^ j ) s = i ; . . . ; m; (2.9) 

In dimension D (with gravitational constant k 2 ), the set of Einstein equations 
R^-RS^ 1 2 = K 2 T$f can be written as R% = K 2 {T$f -T5^ / (D -2)]. Furthermore, like the 
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multidimensional geometry itself, these equations decompose blockwise to Rq — R/2 = k 2 T ° 
and R™* = k 2 [T™* — T5^/(D — 2)]. Using the previous formulas, we obtain 

1 n 

- J2 G ij x t x J + V = 0, (2.10) 

x exp $W(rr)-2( 7 -7o) . (2.11) 

Here, 

Gjj — diSij — d^dj (2.12) 
are the components of the minisuperspace metric, 

m 

V = k 2 Y, A {s) exp [& s \x) - 2( 7 - 7o )] . (2.13) 

3=1 

The integrals (2.9) are used to replace the densities in (2.10), (2.11) by the functions 

It is not difficult to verify that after the gauge fixing 7 = F(x 1 , . . . ,x n ) the Einstein 
equations (2.11) may be considered as the Lagrange- Euler equations obtained from the La- 
grangian 

L = e^-T (± p Gijx'iP - (2.14) 

under the zero-energy constraint (2.10). 

Further we develope an integration procedure which is based on the n-dimensional 
Minkowsky-like geometry. Let IR™ be the real vector space and d, . . . ,e n be the canoni- 
cal basis in H n (i.e. e x = (1, 0, . . . , 0) etc). Let us define a symmetrical bilinear form (, ) on 
]R n by 

(e,, ej) = Sijdj - didj = G^. (2.15) 
The form is nondegenerate and inverse matrix to (Gy) has components 

5 ij 1 

g " = 1: + 2^d- < 2 ' 16 > 

This form < ., . > endows the space IR n with a metric which signature is (—,+,...,+) [19]. 
With this in mind, a vector y e IR n is timelike, spacelike or isotropic, if (y, y) takes negative, 
positive or null values respectively and two vectors y and z are orthogonal if (y, z) = 0. 
Hereafter, we use the following vectors 

x = x 1 (t)e 1 + ... + x n (t)e n , (2.17) 
u = t^d + . . . + u n e n , u % = — J 1 , Ui = d h (2.18) 
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where covariant coordinates Ui of the vector u are introduced by the usual way. Moreover, 
we obtain 

n n 

]T Gijx'x* = (x, x) = di(^) 2 - To. 

i,j=l i=l 

(u,x) = 7o , (u,u) = -753^- 

3 Exact solutions 

Now we suppose that all components are isotropic fluids, i.e. pressures . . . ,pffl in the 
factor spaces M 1? . . . , M n are equivalent for each fluid component. From the mathematical 
point of view it means that 

& s) (x) =fW ( 7o ), s = l,...,m, (3.1) 

where is an arbitrary smooth function on 1R (70 is defined by (2.3)) and variable 
barotropic parameters obtained from (2.7) have the form h\ s \x) = dF^(7 )/d7o for all 
% — 1 . . . , n. 

We use the orthogonal basis 

,f 2 ,...,f n eIR/\ (3.2) 



(u, u) 

where the vector u was introduced by equation (2.18). The orthogonality property reads 

(u, fj) = 0, (fj, f fc ) = S jk , (j,k = 2,...,n). (3.3) 

Let us note that basis vectors f 2 , . . . , f n are space-like, since they are orthogonal to the 
time-like vector u. The vector x e H n decomposes as follows 

u n 

x = 7o 7^\ + ^ (x ' f ^ )f ^ (3 ' 4) 

\u, u; j =2 

We need to obtain the coordinates 70, (x, fi), ... , (x, f n ) as functions of the time t. 
Using the time gauge of the type 

7 = /(To), (3.5) 
we come to the Lagrangian (2.14) in the terms of the coordinates in such basis as follows 

1 / • 2 n \ m 

L = V~^°) + J2(±, f,) 2 -k 2 J2 A^e^-^ +F(s) ^\ (3.6) 

2 y (u, u) - =2 J s=1 

The Lagrange- Euler equations for the coordinates (x, fj), ... , (x, f n ) 

d 



dt 



e 7o-/(7o)^f^j = (3.7) 

5 



give immediatly the following integrals of motion 

e»-/(-»)(x,f,) = a'', j = 2,...,n, 



(3.8) 



where a J is an arbitrary constant. Using (3.8) we may present the zero-energy constraint in 
the form 



/ n m \ 

7 2 + (u, u >e 2 ^^)-^ J2 { ai ) + 2 ^ 2 E A^e F(s) ^) = . 



(3.9) 



s=l 



The last equation admits obtaining of the unknown function 70 in quadratures for arbitrary 
functions /(70), F^^o), . . . , F^(7 ), then the model is integrable. 

To present exact solutions in the Kasner-like form we introduce the following vector 



3=2 



(3.10) 



i=i 



Owing to the orthogonality property given by equation (3.3) the coordinates s l of this vector 
in the canonical basis e 1; . . . , e n satisfy the following constraints 



(s, s) = ]T G y sv = E ( a3 ) 2 - u > = E = 

i,j=l j=2 



(3.11) 



i=l 



Then the vector x we need to find may be presented as 



u 



x = 70 



(u, u) 



+ 



* gIlfa) / ( (s, s) + 2. 2 ± A^e^A d 7o 

y/-(u,u) J V , i / 



s. 



(3.12) 



Here we remind the multi-dimensional generalization of the well-known Kasner solution 
[18]. It reads (for the synchronous time t s ) as follows 



(3.13) 



i=i 



Such a metric describes the evolution of the vacuum model under consideration. Kasner 
parameters e l satisfy the constraints 



11 11 2 



(3.14) 



i=i 



One may easily verify that the introduced parameters and the Kasner parameters e % are 
connected by 



, i = l,...,n. 



Finally, we present the logarithms of scale factors in the Kasner-like form 

-1/2 



(3.15) 



x' =70 



u,u) v 7 v s ti / V ( u > u ) 



(3.16) 
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We recall that coordinates of the vector u are defined by (2.18). The function 70 (t) is 
determined by (3.9). If we use it as the time coordinate the exact solution for the metric 
(2.1) looks as follows 

D — 2 e 270 n 
ds 2 = -— n-ETT^ • d7n + Y exp[2x l ]ds 2 . (3.17) 

4 Examples 

4.1 4-dimensional model 

Now we consider a special case with D = 4 manifold 

M 4 = E x Mj x M\ x M\ 

and a single perfect fluid. The dominant energy condition, applied to the energy-momentum 
tensor (2.5) reduces to the following inequalities for the fluid variable barotropic parameter 
h = F\ lQ ) 

< F'( 7o ) < 2. 
The exact solution for this special model looks as follows 

x* = I + J (<s, s) + 2K 2 Ae F ^)~ 1/2 d 7o • (V - , (4.1) 

p(t) = Aexp [F( 7o (0) - 27o(x)] , p(0 = [1 - F'( 7o (t))] p(t), (4.2) 

where the function 70 (t) = x 1 (t) + x 2 (t) + x 3 (t) is the solution to the separable ordinary 
differential equation 

exp[7o-/(7o)]d 7 o _ [2 dt ^ 



>, s) + 2K 2 Ae F ^) " V 3 

We recall that function /(70) defines a time gauge (/(70) = corresponds to the synchronous 
time). The Kasner parameters e % obey the relations (3.14) with d\ = (I2 = d 3 = 1. The 
nonnegative parameter (s, s) is arbitrary. When (s, s) = one gets an isotropic solution, 
otherwise the solution describes an anisotropic behaviour. 

Let us analyze the behaviour of the exact solution assuming that the function F(j Q ) is 
monotonically increasing with F'(7o) > 0. We choose the following time gauge 

e 70-/(70) = j p'( 70 ) e ^(70) > Q 

Then the equation (4.3) immediately gives 



e^) = i^(t-t ). (t-t 0+ V^|. (4 .4) 
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It follows from (4.4) that F("y Q ) — > +00 and 70 — > +00 as i — > +00. Then the equation (4.1) 
presented in the form 

^UM , , , (4.5) 

d 7o 3 /( s> s ) + 2 K 2 Ae F M 



shows an isotropic expansion of the universe with x % = 7 /3+const, i — 1,2, 3, at late times. 

At earlier times as t — > to one gets e F ^ 70 ^ — > +0 from the equation (4.4). Then the 
equations (4. 3), (4. 5) lead to exp[x l ] ~ (t s — t s o) £l in the main order as the synchronous time 
t s t s o + 0. This shows the Kasner-like behavior near the initial singularity. 

Hereafter we consider this single perfect fluid with the variable barotropic equation of 
state as an associated description of two different perfect fluids: normal matter with positive 
energy density p m and nonnegative pressure p m = (1 — h m )p m and quintessence with positive 
energy density pq and negative pressure pq = (1 — Hq)pq. The description implies 

P = pm + PQ, P = Pm+PQ- (4.6) 

If barotropic parameters h m , hq and h = F'(j ) are specified, we get from the equations 
(4.2),(4.6) 

Pm = A ■ kQ ~ F } l0) ■ exp [F( 7o ) - 2 7o ] , 
riQ — rim 

F'(7o)-/i m lT?( s , 

Pq = A ■ — exp [F( 7o ) - 2 7o ] . 

iiq — ti m 

Then, the energy densities ratio looks as 

Pq _ F'fro) - h m 



a 



pm h Q - F'(7q) 



In what follows we suppose that the barotropic parameters h m and Hq are constant and such 
that 

0<h m <l, Kh Q <2. 
Then the barotropic parameter of the fluid 

h = F'( 7 o) = h Q - 

a(7o) + 1 

is positive and the energy densities ratio a is a function of 70. As we have already shown 
such model describes an isotropic expansion of the universe at late times. If we suppose 
that the model is asymptotically coherent, i.e. the energy densities ratio is asymptotically 
constant 

lim a (70) = olq = const, 

70— >+oo 
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the scale factor and the density are in the main order 



2 _ 

rpi 3 I" '"^ ' CK Q + 1 / _ O 

e ~ t s , p ~ t s as t s — ► +oo, 

where t s is the synchronous time. This is the late times isotropic asymptotic of the solution. 
Further, if the following condition holds 

h-Q — h m 4 

r~r q "q? 

a + 1 3 

the exact solution describes the power law isotropic accelerated expansion at late times. If 
the normal matter is dust (h m = 1) and a according to the observations is equal to 7/3, 
then we obtain Kq > 31/21 or p Q < —10/21 • pq. 

4.2 Multidimensional model 

Now we study the behaviour of the exact solution for the manifold 

M = H x Ml x M 2 d2 x ... x M% n , 

of the dimension D = 4 + d, where Mf is the external flat 3-dimensional space and 
M2 2 , • • • , M^ n are internal spaces. Let us admit that on some early stage of the evolu- 
tion the functions ^^(70), s — 1, . . . ,m, determining the barotropic parameters h\ s \x) = 
dF( s )(7 )/d7o, s = 1, . . . , m and % — 1, . . . , n, of the fluid components, obey the relation 

m A 

A ^e F(a)( ^ = — 2 = const > 0. (4.7) 

s=l 2/t 

In terms of densities p( s \ s — 1, . . . , m, the formula (4.7) looks as follows 

From the physical point of view this relation means that the multicomponent perfect fluid 
appears as a whole as the stiff matter. In this case the exact solution for the scale factors is 

= al-(t s -t s0 y\ (4.8) 

where 



D-l V D-lJ \ (8,8) J 

Using the Kasner-like constraints (3.14) one easily obtains the similar constraints for the 
parameters e l 

n n „ 1 p< 9/ 4 \ — 1 



i=l i=l 
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It is readily seen from (4.9) that if the parameter A/(s,s) is small, then 5 m 1 and, conse- 
quently, e l ~ e l . Then the model describes the Kasner-like behavior of the multidimensional 
model with the expansion of one part of the spaces Mf, M$ 2 , . . . , M^ n and the contraction 
of another part. From the physical viewpoint the behavior with expansion of the external 
space Mf and the contraction of internal space (or spaces) is of the most interest. Let us 
analyze this solution for the space-time manifold M 4+d = Rx Mf x M|, i.e. when there is 
only one <i-dimensional internal space M|. In this case the solution to the set of equations 
(4.9) reads 



1 ± 



d 



l + 2/d 

d + 3 3 Y 1 + ^/(s,s 



3- 



e 2 = 



it 



3- 



l + 2/d 



d + 3\ ^\ l + A/(s,s)' 



(4.10) 



We notice that within the model the function 



JJBD) 



const • (t s —t 



is the scale factor of the external space with respect to the Jordan-Brance-Dicke frame (see, 
for instance, [22], [23], [24]). With respect to the Einstein frame the scale factor of the external 
space reads 



(E) _ (JBD) d . x 2/ 2 



= const • (t s - t s0 ) £ , 



where 



e 1 = e 1 + d ■ e 2 /2 



d + 2 d 



d+3\ 2 6 \ l + A/(s,s 



3- 



1 + 2/d 



< e 1 



(4.11) 



In order to obtain the contraction of the internal space (i.e. e 2 < 0) we take the upper signs 
in the formulas (4. 10), (4. 11) with the following obvious condition 



1 + 2/d 



> 1. 



\ l + A/(s,s) 

Then A/ (s, s) e (0,2 + 6/d). Using the last one easily gets the following inequalities 



1 -i 1 + 1 
- < e 1 < - 



3(l + 2/d) d + 2-l-j3(l + 2/d) l 

< e < -. 



d + 3 



d + 3 



The analysis shows that if the solution describes contraction of the internal space then the 
external space in this solution may expand only with deceleration. 
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